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This is an addendum to the previous publication, P. Ko and J.-h. Park, Phys. Rev. D80, 035019 
(2009). The semileptonic charge asymmetry in B s decays is discussed in the context of general 
MSSM with gluino-mediated flavor and CP violation in light of the recent measurements at the 
Tevatron. 
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In this addendum to Ref. [1 j, we discuss the semilep- 
tonic charge asymmetry in the B s decays in general 
SUSY models with gluino-mediated flavor and CP vio- 
lation, in light of the recent measurements of like-sign 
dimuon charge asymmetry by D0 Collaboration at the 
Tevatron. The model is described in Ref. [1 , to which we 
refer for the details of the model and other phenomeno- 
logical aspects related with B s —B s mixing, the branching 
ratio of and CP asymmetry in B — >• X s 7, Bd —> 4>Ks and 
CP asymmetry in B s — >> J/ijxj). 

One can define the semileptonic charge asymmetry in 
the decay of B q mesons as 



in the asymmetry, 



r(B$(t) -> ii+ x) + r(B°(t) -+ ti-x) ' 



(i) 



for q = s. In terms of the matrix elements of the 
effective Hamiltonian describing the damped oscillation 
between B® and B® the asymmetry a q sl is given by 
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where (j) q = arg(— Mj^/r^). That is, this is another ob- 
servable measuring CP violation in B q -B q mixing. We 
take the approximation, T\ 2 = T^ 2 SM , since the leading 
contribution comes from the absorptive part of the box 
diagrams for B q -B q mixing and there is no new com- 
mon final state into which both B q and B q can decay 
in our scenario. The size of M\ 2 is fixed by the AM q 
data up to hadronic uncertainties. Then, can be re- 
garded as a sine function of qb qi multiplied by the factor 
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|Mf 2 |- This curve is traversed as one allows for 



arbitrary supersymmetric contributions to M\ 2 obeying 
the AM q constraint. Combining the SM predictions [2], 
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(49.7 ±9.4) x 10~ 4 , 
(4.2 ±1.4) x 10" 3 , 
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one finds the vanishingly small asymmetry a sl 

io- 5 . 

Recently, the D0 collaboration reported a measure- 
ment of like-sign dimuon charge asymmetry [3]. They 
interpreted the result as coming from the mixing of neu- 
tral B mesons and have found an evidence for an anomaly 
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where + and N b 



are the number of events where de- 
cays of two b hadrons yield two positive and two negative 
muons, respectively. Their result shows a discrepancy of 
3.2<j from the SM expectation. This asymmetry consists 
of coming from B^ decays as well as from B s . One 
can extract the asymmetry relevant to the B s meson us- 
ing the measured value of and the result by D0 is 



-0.0146 ±0.0075. 



(5) 



This is 1.9<j away from the SM prediction. We shall use 
this data in the following discussion. 

This D0 result has drawn interest in new physics ex- 
planations [1H8]. (For earlier works, see e.g. Refs. [9]- 
HT].) Some of the works consider extra contributions to 
r^ 2 since the dimuon charge asymmetry depends on it as 
well as on M q 2 [5] |6] . This approach also has a possibil- 
ity of altering |Ar s | even though its current experimental 
value is in agreement with the SM one, 2 |T^ 2 SM cos^ M | 
[2j [12j [13]. As we said, r^ 2 is fixed in the present work 
and we are left only with the option of modifying Mf 2 . 
Therefore, |AT S | shall become smaller than its SM pre- 
diction as | (j) s | grows up to 0(1). 

We perform the numerical analysis in the same way as 
in the main article pQ. The crucial ingredient for eval- 
uating is the range of (j) s to be used. Following the 
latest reports from D0 [3] and CDF [14] . there have been 
a couple of attempts to make a global fit of B s -B s mix- 
ing parameters including (j) s [4j[6|. However, the official 
combination is not available yet. Partly because of this 
reason and partly for the sake of coherent presentation, 
we keep using the range used in Refs. [T| [T5], 

(j) s e [-1.10, -0.36] U [-2.77, -2.07]. (6) 



As a matter of fact, this range is not very different from 
the 2a interval found in Ref. [6 . As for r^/M^ M , 
we take its central value from Eqs. (J3|. Considering the 
error in this ratio could add 20% more of uncertainty to 
the thickness of the band in the following figures. 

We show as a function of <p s for tan/3 = 3 in Figs. HI 
The four plots are for the LL, the RR, the LL = RR, and 



2 




(c) LL = RR (d) LL = -RR 

FIG. 1. Plots of ali as a function of <f> 8 for the four different cases with tan/3 = 3. The hatched gray region leads to the lightest 
squark mass < 100 GeV. The hatched region is excluded by the B — » X s j constraint. The light gray region (cyan online) is 
allowed by AM S . The dark gray region (blue online) is allowed both by AM S and <fi s . The black square is the SM point. The 
dashed and solid lines (both red online) mark the la and 2a ranges of aj, respectively. 



the LL = —RR cases, respectively. One can immediately 
notice the aforementioned sinusoidal dependence of on 
<f> 3 , coming from Eq. ^ and the AM S constraint. This 
feature is not only true of all the cases shown here but also 
of any new physics model that does not affect Tf 2 - The 
nonzero thickness of the band arises from the uncertainty 
in AM S . The difference between and its central value 
is at least about 1.0a. This discrepancy becomes worse 
but only slightly after <j) s is restricted inside its preferred 
ranges (colored in blue). If one incorporates the B — >• 
X s j constraint, substantial part of the blue regions is 
excluded, in particular in the upper two cases with one 
insertion. Even then, however, the lowest possible value 
of ~ —0.006 within the blue region does not change. 
In the lower two cases with two insertions, B — ^ X s j 
does not play an important role since the supersymmetric 
effect on B s -B s mixing is enhanced. 



Plots for tan/3 = 10 are displayed in Figs. [2] The 
model-independent characteristics dictated by Eq. Q re- 
main exactly the same as in the previous set of figures. 
The only difference is the stronger B — >> X s j constraint 
due to higher tan/3. Here, it excludes more part of the 
blue regions. Again, this is particularly true of the up- 
per two cases in which is restricted closer to its SM 
value. In Fig. |2|a), AM S , </> s , and B — >> X.7, together 
allow to be as low as —0.003. In Fig. |2^b), there is no 
solution satisfying all the three constraints. One could 
get ~ —0.0006 if <p s were not limited. In the lower 
two cases, the lowest a| 1? compatible with AM S and </> s , 
is almost the same as in Figs, [l] 

We summarize. We have examined how is in- 
fluenced by the LL and/or RR mass insertions. For 
tan/3 = 3, one can reduce the discrepancy between a s sl 
and its SM expectation from 1.9a down to 1.0a in each 
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FIG. 2. Plots with tan/3 = 10. The meai 

of the LL, RR, LL = RR, and LL = —RR cases, obeying 
the AM S , B — » X s 7, and (j) s constraints. This amounts 
to reduction of the A b sl tension from 3.2<j down to 2.2<j 
if one assumes no new physics in the b — >• d transition. 
For tan/3 = 10, it becomes difficult for the LL and RR 
cases whereas the LL = and LL = — iU? cases are 
less limited by B — >• X s 7. 
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(b) RR 
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(d) LL = -RR 
g of each region is the same as in Figs, [I] 

NOTE ADDED 

While we were waiting for the approval for submission, 
a paper by J. K. Parry appeared on the e-print archive 
that employs a related model [8 . However, the flavor 
structure of the squark mass matrix therein is different 
from any of those here. As far as squarks are concerned, 
he considers only one case where (^3)^^ is a variable 
parameter and (S^) ll is fixed to a value that comes from 
renormalization group running. This way of parameter 
scan is not covered in this work. He does not display the 
B — >• X s 7 constraint on his plots, but it may not be very 
restrictive in his case depending on \i and tan/3. (See e.g. 
Fig. 4 in Ref. [16].) 



[1] P. Ko and J.-h. Park, Phys. Rev. D 80, 035019 (2009) [2] A. Lenz and U. Nie rste, JHEP 0706, 072 (2007) 
|arXiv:0809.0"705| [hep-ph] ] . |arXiv:hep-ph/0612167] . 



4 



[3] V. M. Abazov et al [DO Collaboration], Phys. Rev. D 
82, 032001 (2010) [arXiv: 1005.2757] [hep-ex]]. 

[4] Z. Ligeti, M. Papucci, G. Perez and J. Zupa n, Phys. Rev. 
Lett. 105, 131601 (2010) arXiv: 100 6.0432 [hep-ph]]. 

[5] A. Dighe, A. Kund u and S. N andi, Phys. Rev. D 82, 
031502 (2010) [arXiv: 1005.4051] [hep-ph]]. 

[6] C. W. Bauer and N. D. Dunn, |arXiv: 1 006. 1629 [hep-ph]. 

[7] B. A. Dobrescu, P. J. Fox a nd A. Martin, P hys. Rev. 
Lett. 105, 041801 (2010) |arXiv:1005.4~238l [hep-ph]]; 
A. J. Buras, M. V. Carlucci, S. Gori and G. Isidori, JHEP 
1010, 009 (201 0) |arXiv:1005.5310| [hep-ph]1; C. H. Chen 
and G. Faisel, arXiv:1005.4582| [he p-ph]; D. Choudhury 
and D. K. Ghosh, |arXiv:1006.2171| [hep-ph] ; N. G. Desh- 
pande, X. G. He and G. V alencia, Phys. Rev. D 82, 
056013 (2010) [arXiv: 1006. 1682] [hep-ph]]; C. H. Chen, 
C. Q. Geng and W. Wang, JHEP 1011, 089 (2010) 
[arXiv: 1006.52161 [hep-ph]]. 

[8] J. K. Parry, |arXiv: 1006.5331] [hep-ph] . 

[9] G. L. Kane, P. Ko, H. b. Wang, C. Kolda, J.-h. Park 
and L. T. Wang, Phys. Rev. Lett. 90, 141803 (2003) 



arXiv:hep-ph/0304239 ; Phys. Rev. D 70, 035015 (2004) 

arXiv:hep-ph/02 1 2092] . 
[10] A. Lenz, Phy s. Rev. D 76, 065006 (2007) 

|arXiv:0707.1535| [hep-ph]] . 
[11] L. Randall and S. f. S u, Nucl. Phys. B 540, 37 (1999) 

arXiv:hep-ph/9807377 ; K. Kawashima, J. Kubo and 
A. Lenz, Phys. Lett. B 681, 60 (2009) |arXiv:0907.2302| 

[hep-ph]]. 

[12] E. Barberio et al [Heavy Flavor Averaging Group], 

|arXiv:0808.1297| [hep-ex] . 
[13] S. Esen et al, P hys. Rev. Lett. 105, 201802 (2010) 

|arXiv:1005.5T77| [hep-ex]] . 
[14] L. Oakes (CDF Coll aboration), talk at FPCP 2010, May 

25-29, Torino, Italy, http://agenda.infn.it/getFile. 

py/access?contribld=12&resld=0&materialld= 

slides&conf Id=2635 
[15] M. Bona et al [UTfit Collaboration], PMC Phys. A 3, 6 

(2009) |arXiv:0803.0"659| [hep-ph]]. 
[16] P. Ko, J.-h. Park and M . Yamaguchi, JHEP 0811, 051 

(2008) arXiv:0809.2784 [hep-ph]]. 



